AD-A268  512 

liiiiilflii 


UNIVERSITY 
OF  WISCONSIN 


CENTER  FOR  THE 
MATHEMATICAL 
SCIENCES 


Center  for  the  Mathematical  Sciences 
University  of  Wisconsin — Madison 
1308  W.  Dayton  Street 
Madison,  Wisconsin  53715-1149 


DtlC 

ELECTE 

ftUGX7  1993 

E 


D 


July  1993 


(Received  July  27,  1993) 


Approved  for  public  release 
Distribution  unlimited 


Sponsored  by 

National  Science  Foundation  Office  of  Naval  Research 

1800  G  Street,  N.W.  Arlington,  VA  22217 

Washington,  DC  20550 


8  '16  0 


UNIVERSITY  OF  WISCONSIN-MADISON 
CENTER  FOR  THE  MATHEMATICAL  SCIENCES 


WAVE  STRUCTURE  INDUCED  BY  FLUID  DYNAMIC  LIMITS 
IN  THE  BROADWELL  MODEL 

Athanasios  E.  Tzavaras 
CMS  Technical  Summary  Report  ^94-2 
July  1993 

Abstract.  Consider  the  fluid  dynamic  limit  problem  for  the  Broadwell  system  of 
the  kinetic  theory  of  gases,  for  Riemannn,  Maxwellian  initial  data.  The  formal  limit  is 
the  Riemann  problem  for  a  pair  of  conservatioir  laws  and  is  invariant  under  dilations 
of  coordinates.  The  approach  of  self-similar  fluid  dynamic  limits  consists  in  replacing 
the  mean  free  path  in  the  Broadwell  model  so  that  the  resulting  problem  preserves  the 
invariance  under  dilations.  The  limiting  procedure  was  justified  in  [ST].  Here,  we  study  the 
structure  of  the  emerging  solutions.  We  show  that  they  consist  of  two  wave  fans  separated 
by  a  constant  state.  Each  wave  fan  is  associated  with  one  of  the  characteristic  fields  and 
is  either  a  rarefaction  wave  or  a  shock  wave.  The  shocks  satisfy  the  Lax  shock  conditions 
and  have  the  internal  structure  of  a  Broadwell  shock  profile. 


AMS  (MOS)  Subject  Classifications;  35L60,  35Q20,  76P05 

Key  Words  and  Phrases:  Broadwell  model,  fluid  dynamic  limit,  shock  waves^ 

I  Justification 


I 


OTIC  QUALITY  mSPECTED  3 


By _ 

Distribution/ 


Availability  Codes 

Dist 

6± 

Avail 

Spe 

and /or 
cial 

Research  partially  supported  by  the  National  Science  Foundation  under  Grant  DMS- 
9209049  and  the  Office  of  Naval  Research  under  contract  N00014-93-0015. 


WAVE  STRUCTURE  INDUCED  BY  FLUID  DYNAMIC  LIMITS 
IN  THE  BROADWELL  MODEL 


§1.  Introduction 


Being  among  the  simplest  models  in  the  kinetic  theory  of  gases,  the  Broadwell  model  has 
served  as  a  paradigm  to  understand  the  phenomenon  of  relaxation  and  the  transition  from 
a  microscopic  to  a  macroscopic  description  of  gases.  It  consists  of  the  system  of  seinilinear 
hyperbolic  equations 


dt 


(1.1) 


dt 


and  is  associated  with  a  six-velocity  model  describing  a  system  of  particles  with  identical 
masses  that  move  along  three  mutually  orthogonal  directions  with  speeds  ±c  and  interact 
according  to  a  mechanism  of  equiprobable  binary  collisions  (Broadwell  [B]).  The  system 
(1.1)  derives  from  the  six-velocity  model  when  specializing  to  one- dimensional  flows,  for 
which  the  densities  of  particles  moving  in  the  directions  orthogonal  to  the  flow  are  all 
equal.  (We  refer  to  Platkowski  and  Illner  [PI]  for  the  derivation  and  references). 

In  this  context  the  function  /  =  (/i,/2,/3),  defined  for  (x,<)  G  IR  x  IR'*',  describes  the 
densities  of  particles:  f\  for  particles  moving  in  the  positive  z-direction,  /2  in  the  negative 
z-direction  and  /s  in  each  of  the  positive  or  negative  y-  or  z-directions.  The  speed  c  is 
taken  here  equal  to  one  and  £  stands  for  the  mean  free  path,  a  measure  of  the  average 
distance  between  successive  collisions.  The  right  hand  side  is  called  the  collision  operator 


and  measures  the  rate  of  gain  (or  loss)  in  densities  of  particles  effected  through  collisions. 
It  is  characterized  by  the  quantity 


W)  =  /|-/i/2.  (1-2) 

The  zeroes  of  Q{f)  axe  the  states  of  equilibrium  for  the  system,  /I  =  /1/2,  and  are  called 
Maxwellians.  Finally,  associated  with  each  /  are  the  quantities 

Pf  =  f\  +  f2  +  4/3  ,  mf  =  fi  —  f2  (1.3) 

measuring  the  local  density  and  momentum  flux  in  the  x— direction,  respectively. 

The  limit  when  the  mean  free  path  approaches  zero  is  known  as  the  fluid  dynamic 
limit.  For  small  mean  free  path  the  strong  interactions  of  particles  entail  a  macroscojjic 
description  of  the  flow  to  become  meaningful.  In  the  case  of  the  Broadwell  model  the 
induced  macroscopic  ’’Euler  equations”  are  easy  to  identify.  First  rewrite  (1.1)  in  the  form 
^(/i  4-/2+  4/3)  +  ^(/i  -  /z)  =  0 , 

+  =  (1.4) 

Formally,  eis  e  — >  0,  it  is  expected  that  the  first  two  equations  pass  in  the  limit,  while  the 
third  causes  the  limiting  /  to  be  a  local  Maxwellian.  Thus  the  limit  fluid  equations  become 

|(/l  +  /2  +4(/,/2)''")  4-  ^(/.  -  /2)  =  0, 

^(/i  -/2)+^(/l  +/2)=0. 

The  corresponding  macroscopic  density  and  momentum  of  the  fluid  are  given  by 

P  = /i +4(/i/2)'^^ +/2,  m  =  pu  =  /i-/2.  (1.6) 
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The  algebraic  system  can  be  easily  inverted  and  leads  to  an  alternative  form  of  the  limit 
’’Euler  equations”,  in  terms  of  the  macroscopic  variables  (p,  u), 


dt 


+  gj(pu)  =  0, 


|(P«)  +  ^(P!7(«))  =  0, 


(1.7) 


where  g{u)  :=  +  Zu^yl'^  —  1  .  They  form  a  strictly  hyperbolic,  genuinely  nonlinear 

system  of  conservation  laws  (Caflisch  [C]). 

The  justification  of  the  fluid  dynamic  limit  has  been  the  object  of  several  investigations. 
We  refer  to  Cercignani  [Ce]  for  a  survey  of  the  literature  on  the  Boltzmann  equation 
and  to  Platkowski  and  Illner  [PI]  for  result?  on  discrete  velocity  models  of  kinetic  theory. 
For  the  Broadwell  model,  the  fluid  dynamic  limit  is  understood  for  smooth  solutions  of 
the  limit  ’’Euler  equations”  (Inoue  and  Nishida  [NI],  Caflisch  and  Papanicolaou  [CP]). 
Regarding  the  case  of  solutions  with  shocks,  we  mention  the  studies  on  stability  (in  time)  of 
traveling  wave  solutions  (Kawashima  and  Matsuniura  [KM])  or  rarefaction  wave  solutions 
(Matsumma  [M])  for  the  Broadwell  model,  and  a  recent  study  by  Xin  [X],  showing  that 
a  given  piecewise  smooth  solution  with  noninteracting  shocks  of  the  limit  fluid  equations 
can  be  approximated  by  solutions  of  the  Broadwell  system  as  £  -»  0,  that  gives  a  definitive 
answer  to  one  direction  of  the  problem.  The  converse  problem,  to  show  that  a  given 
family  of  solutions  to  the  Broadwell  system  converges  globally  in  time  to  a  fluid-dynamical 
solution,  remains  at  present  open. 

Iii.sight  in  the  latter  direction  is  provided  by  the  approach  of  self-similar  fluid  dynamic 
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limits  (Slemrod  and  Tzavaxas  [ST]).  For  Riemann,  Maxwellian  data 


:.0)  =  I  ^  °  with  =  (ffJt.fi') 


QiD  =  0(/+)  =  0 .  /,* ,  ff ,  A*  >  0 , 


±  r±  f± 


the  solutions  of  the  limit  fluid  equations  (1.5)  are  expected  to  be  self-similar  functions 
of  if  =  x/t.  On  the  other  hand,  the  Broadwell  system  does  not  possess  space-time  dila- 
tional  invariance  and  does  not  admit  self-similar  solutions  of  that  type.  Motivated  by  an 
cuialogous  idea  for  systems  of  conservation  laws  (Dafermos  [Di]),  one  considers  a  modified 


Broadwell  system 


dfi  dfi  _  1  ^ 

dt  dx 
dt 


which  does  preserve  the  dilational  invariance  {x,t)  —*  {ax,  at),  a  >  0.  The  resulting 
Riemann  problem  (1.8  —  1.9)  admits  self-.similar  solutions  of  the  form  /  =  f{x/t).  These 
are  constructed  by  solving  the  singular  boundary  value  problem 

-(e-i)/i'(e)  =  (/l-/i/2)/£ 


-(e+l)/2(0=(/|-/l/2)/£ 

-e/3=-(/|-/l/2)/2£ 


(1.10) 


/(-i)  =  r,  /(+i)  =  r, 

for  f  €  [~1»1]  and  subject  to  (M),  referred  to  as  problem  {Ve).  It  has  the  following 
property  [STj:  Any  family  of  solutions  {/''}e>o  to  (Pe)  corresponding  to  fixed  data  is 
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of  uniformly  bounded  variation  and,  as  a  consequence,  there  exists  a  subsequence  {/^"  } 
with  £„  — ♦  0  and  a  function  of  bounded  v-ariation  /  such  that  /*'■  — ^  /  pointwise.  The 
limiting  /  is  a  local  Maxwellian,  satisfying  /s  =  a.e  in  [—1,1]-  If  /  is  extended 

to  the  entire  real  line,  by  setting  /  =  /“  on  (-00,  -1)  and  /  =  /'*'  on  (1,  oc),  the  function 
f{xjt)  becomes  a  weak  solution  of  the  Riemami  problem  (1.5,  1.8). 

The  goal  of  this  article  is  to  investigate  the  .structure  of  the  limiting  function  /  con¬ 
structed  via  self-similar  fluid  dynamic  limits.  Of  special  interest  is  to  clarify  the  behavior 
of  /(O  at  points  of  discontinuity.  The  Broadwell  model,  being  among  the  simplest  in  ki¬ 
netic  theory,  is  a  case  study  to  understand  the  mechanism  of  relaxation  and  the  associated 
admissibility  restrictions  imposed  on  shocks. 

The  presence  of  relaxation  mechanisms  is  natural  in  many  physical  contexts  and  has  been 
investigated  extensively  both  for  specific  models,  such  as  Broadwell,  but  also  for  systems 
of  two  conservation  laws  with  relaxation  {e.g.  Liu  [L]).  For  (1.1),  the  study  of  shock 
profiles  by  Broadwell  [B]  captures  the  regularizing  effect  induced  on  shocks  by  rela.xation. 
These  shock  profiles  have  been  compared,  and  are  in  close  agreement  for  weak  shocks,  with 
traveling  wave  solutions  of  an  associated  system  of  viscous  conservation  laws  arising  from 
(1.1)  via  the  Chapman-Enskog  expansion  (Caflish  [C]).  An  important  difference  of  the  self¬ 
similar  relaxation  investigated  here  is  that  it  penalizes  the  whole  wave  fan  simultaneously. 
Comparisons  between  the  structure  entailed  by  self-similar  limits  and  the  Broadwell  shock 
profiles  are  carried  out  in  the  text. 

The  structure  of  the  limiting  solution  is  completely  characterized  via  this  approach.  It 


turns  out  that  /  consists  of  two  wave  fans  separated  by  a  constant  state.  Each  wave  fan  is 
associated  with  one  of  the  characteristic  fields  of  the  limit  fluid  equations  and  consists  of 
either  a  single  rarefaction  or  a  single  shock.  The  shocks  satisfy  the  Lax  shock  conditions 
and  have  the  internal  structure  of  a  Broadwell  shock  profile. 

The  approach  of  self-similar  viscous  limits  to  study  the  structure  of  the  Riemann  problem 
was  initiated  by  Dafermos  [D2],  in  the  context  of  systems  of  two  conservation  laws.  Several 
of  the  ideas  developed  in  [Do]  are  used  here,  but  the  basic  strategy  is  different.  One  starts 
with  some  representation  formulas  expressing  the  rate  of  colli.sions  Q{f^)/e  (and  thus  also 
the  derivatives  fj')  as  an  averaging  process  depending  on  functions  of  the  solution  f^. 
Using  the  estimates  on  the  total  variation  of  from  [ST],  it  can  be  shown  that  there  is  a 
finite  Borel-Stieltjes  measure  so  that  along  a  subsequence 

<?(/*) 

- ^  u  weak-*  in  measures . 

£ 

The  measure  u  incorporates  the  form  of  the  limiting  /.  Although  u  can  not  be  explicitly 
computed  via  this  approach,  it  is  possible  to  establish  several  properties  of  u.  Its  support 
coincides  with  the  set  where  /  is  nonconstant.  Most  important,  using  the  representation 
formulas,  it  turns  out  that  at  points  of  suppi/  a  function  g  related  to  the  antiderivatives 
of  the  eigenvalues  is  maximized  (Proposition  4.4).  The  maximization  property  captures  (a 
large  part  of)  the  restrictions  imposed  on  solutions  by  the  relaxation  process.  It  provides 
the  behavior  of  the  eigenvalues  to  the  limit  fluid  equations  on  suppi/  (Corollary  4.5).  The 
first  part  of  the  proof,  carried  out  in  Sections  3  and  4,  is  completed  by  characterizing  the 
behavior  of  right  and  left  derivatives  at  points  in  the  wave  fans  (Proposition  4.6). 
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The  emerging  strvicture  is  that  typical  (at  least  for  small  data)  of  strictly  hyperbolic 
blit  not  necessarily  genuinely  nonlinear  systems.  The  reason  is  that  although  the  analysis 
in  Sections  3  and  4  makes  essential  use  of  the  strict  hyperbolicity  property  of  the  limit 
fluid  equations,  it  never  uses  its  genuine  nonlinearity  properties.  Complementing  this 
picture,  the  geometry  of  the  shock  curves  for  (1.5)  is  employed  in  Section  6  to  rule  out  the 
possibihty  of  contact  discontinuities  and  arrive  at  the  stated  result.  Finally  in  Section  5, 
an  idea  from  [D2]  is  adapted  to  show  that  at  points  of  discontinuity  the  family  has  the 
internal  structure  of  a  Broadwell  shock  profile  (Proposition  5.1). 

§2.  Preliminaries.  The  boundary  value  problem  {Vs)- 

First  we  review  some  preliminary  information  on  solutions  of  the  boundary  value  problem 
(Vs).  The  material  presented  in  this  section  is  taken  out  of  [ST]. 

As  the  underlying  differential  equations  are  singular,  the  meaning  of  solutions  needs 
clarification.  A  function  /  =  ( /i  ,  /z ,  /a  )  defined  on  [—1,1]  solves  (Vs)  if  fj  €  C[— 1, 1] 
satisfy  the  boundary  conditions  /y(±l)  =  f^,j  =  1,2,3,  and  the  weak  form  of  the 
equations 

-(f-  +  f  Hr)dT  =  i  f  Q{f(T))dr,  (2.1) 

J  a  ^  J  a 

-(«  +  l)/2(«ll  +  f  h(r)dr  =  i  /  (?(/(r))* ,  (2.2) 

J  a  ^  J  a 

+  j'  h(r)dr  =  QU{r))dr ,  (2.3) 

for  all  a,  6  €  [—1,1].  It  follows  immediately  that  a  solution  in  the  above  sense  is  contin¬ 
uously  differentiable  (except  possibly  at  the  singular  points  ±1  and  0)  and  satisfies  the 
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strong  form  of  equations  (1.10)  on  (  —  1,0)  and  (0, 1).  A  more  detailed  analysis  shows  that 
for  positive,  Maxwellian  data: 

(i)  MO  >  0  for  ^  €  [-1, 1] 

(ii)  W(-l))  =  Q(/(0))  =  Q(/(+l))  =  0 

(hi)  fj  €  1,1]  for  some  0  <  ««  <  1,  with  increasing  as  e  decreases;  moreo'^er, 

there  is  £o(/^)  >  0  such  that  for  e  <  Sq  the  functions  fj  are  Lipshitz  continuous  on  [—1,  Ij. 

(i'’)  Q(f)  does  not  change  sign  on  the  intervals  (  —  1,0)  or  (0,1).  As  a  consequence  the 
shapes  of  the  functions  fj  belong  to  one  of  the  following  categories: 

Ci:  Q(f)  >  0  on  (  —  1,0)  and  (0,1);  /i  is  increasing  on  (—1,1),  /2  is  decreasing  on 

(  —  1,1),  fs  is  decreasing  on  (  —  1,0)  and  increasing  on  (0,1). 

C2'  Q{f)  <  0  on  (—1,0)  while  Q{f)  >  0  on  (0,1);  f\  is  decreasing  on  (—1,0)  and  in¬ 

creasing  on  (0, 1),  /2  is  increasing  on  (—1, 0)  and  decreasing  on  (0, 1),  fz  is  increasing 
on  (-1,1). 

C3:  Q{f)  <  0  on  (—1,0)  and  (0,1);  f\  is  decreasing  on  (—1,1),  /2  is  increasing  on 
(—1,1),  /a  is  increasing  on  (—1,0)  and  decreeising  on  (0,1). 

C4:  Q{f)  >  0  on  (—1,0)  while  Q(/)  <  0  on  (0, 1);  f\  is  increasing  on  (—1,0)  and  decreas¬ 
ing  on  (0, 1),  /2  is  decreasing  on  (—1,0)  and  increasing  on  (0,1),  /a  is  decreasing 
on  (-1,1). 

C5:  Q{f)  =  0  on  (-1,0)  eind/or  Q{f)  =  0  on  (0, 1);  In  this  case  /1,/2,/a  are  constant 
on  the  region  where  Q{f)  =  0,  and  /1,/2,/a  have  the  behavior  indicated  in  Cases  1 
-  4  where  Q{f)  >  0  or  Q{f)  <  0. 


8 


Henceforth,  the  data  are  fixed  subject  to  (M)  and  it  is  assumed  that  for  each  £  >  C 
the  problem  (P«)  admits  a  solution  The  question  of  existence  is  considered  in  [ST]. 
The  difficulty  of  this  problem  is  that  the  underlying  equations  are  singular  and  that  one 
needs  to  guarantee  existence  of  solutions  for  every  value  of  the  parameter  £  >  0.  This  has 
only  been  accomplished  for  a  certain  class  of  boundary  data  (c./.  [ST]).  However,  such 
restrictions  do  not  enter  in  any  other  way  in  the  present  analysis. 

The  functions  are  extended  to  the  whole  real  line  by  setting  /'  =  f~  on  {-oo,  -1) 
and  =  /■*■  on  (1.  oo).  The  restrictions  on  the  shapes  of  solutions  are  the  key  ingredient 
of  the  following  theorem  [ST,  Leruna  1.2,  Theorem  2.1]. 

Theorem  2.1.  Let  {/'}£>o  be  a  family  of  extended  solutions  of  (Vs)  corresponding  to 
data  satisfying  (M).  Then: 

(a)  There  are  positive  constants  rtij  ,  Mj  and  Kj ,  j  =  1,2,3,  depending  on  the  boundary 

data  but  independent  of  e  such  that 

0<m,  </;(0<M,  ,  ee[-l,l]  (2.4) 

<  A',  (2.5) 

(b)  There  exists  a  subsequence  {/'"  }  with  £„  -»  0  and  a  positive,  bounded  function  f  of 

bounded  variation  such  that  — ♦  /  pointwise  on  the  reals.  The  function  f  satisfies 

h  =  \//>  h  for  a.e.  ^  €  [-1,1]  (2.6) 


r/  on  (-00,-1] 

\  /+  on  [l,oo) 
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and  the  balance  of  mass  and  momentum  equations 

+  h + ^(f,  -  h)  =  0, 
-  fl)  +  ^(/>  +  /2)  =  0  , 

in  the  sense  of  distributions  and  in  the  sense  of  measures. 

The  function 


f(j)  =  ^  (-00,00)  X  (0,co), 


(2.8) 


is  a  weak  solution  of  the  Riemann  problem  (1.5, 1.8).  Indeed,  (1.8)  is  certainly  satisfied. 
That  the  weak  form  of  (1.5)  is  satisfied  follows  from  the  weak  form  of  (2.7)  by  using  a 
change  of  test  functions. 

§3.  Structure  of  solutions  I 

Let  (/'},>«  be  a  family  of  solutions  to  {Ve)  corresponding  to  fixed  boundary  data 
subject  to  (M).  The  functions  /'  satisfy 


with 


_  _±_  Qun 

i-e  e 

_ ]_  Qif) 

^  +  1  £ 

1  (?(/') 


(3.1) 

(3.2) 

(3.3) 


Q{n=  Qor  =  f^-nn- 


(3.4) 
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The  form  of  the  above  equations  suggests  to  monitor  the  quantity  Q(f^)f£  and  study  its 
limiting  behavior  as  £  — >  0. 


Representation  formulas 

First  certain  representation  formulas  for  Q{f‘)f£  and  the  derivatives  /J'  are  derived. 
We  introduce  the  notation 


,  n  fi  ,  n 

^  =^-r37+T 

and  remark  that  the  function  depends  implicitly  on  e  through  the  dependence  on  the 
solution  /*.  A  simple  computation,  using  (3.1  —  3.3),  shows  that 


dQin  _  1 .  /f 

di  £  ^  +  1 


JL 

1-^ 


+ 


jc'o(r). 


(3.6) 


After  an  integration  (3.6)  leads  to 


Q(/nO) 

£ 


■  !«/'(<.-))  exp{i/'_c'W*) 
;Q(/‘(«+))  exp  I 


for  ^  G  (-1,0) 
for  ^  G  (0, 1) 


(3.7) 


where  a_,  0+  are  any  two  points  with  —  1  <  a-  <  0  <  «+  <1.  It  is  convenient  to 


introduce  a  notation 


s‘«)  =  < 


C‘(s)ds 


«  €(-1,0) 

?  €(0,1) 


(3.S) 


for  the  integrals  appearing  in  (3.7). 
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Integrating  (3.1)  and  using  (3.7),  we  obtain 

Q(/'(«-)) 


fm  -  = 

/.■*■  -  sm  = 


£  ^0  1  ~  C 


(3.9) 


In  turn,  substitution  of  (3.7)  and  (3.9)  in  (3.1)  gives 


/f'(0  =  (/f(o)-A-) 

fi'io  =  {ft  -  fm) 


f]  dC  ’ 


ee(-i,o), 
e  e(0,i). 


(3.10) 

(3.11) 


Relations  (3.10  —  3.11)  are  conceptually  interesting  as  they  express  /f'  as  an  averaging 


jarocess.  Note  that  the  weighting  factors  depend  on  the  solution  /' .  Analogous  relations 
hold  for  the  derivatives  and  f^'. 

In  a  similar  manner  using  (3.1  —  3.3)  and  (3.7),  we  obtain 


<?(/*(<>-)) 


fm  -  fr 


fm-ff 


fm  -  /r 


<?(/•(«+»  /,+  -/,'(o)  _  - /?(o)  ft -fm 


(3.12) 

(3.13) 


^ 

Equations  (3.12  —  3.13),  in  conjunction  with  (3.7)  and  (3.1  —  3.3),  show  that  Q{f^)/e  and 
the  derivatives  may  be  expressed  in  a  variety  of  ways.  In  the  sequel  we  work  with 


Qin  _  f  (/f(o)-/r)Mi 

=  I  (/«0) -/,*)#<?. 

where  n‘_  is  defined  on  the  interval  (  —  1,0)  by 


= 


ei9‘(0 


on  (  —  1,0) 
on  (0, 1) 


J  11 


(3.14) 


J-i  i-C 


(3.15) 


dC 
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)u^  is  defined  on  the  interval  (0, 1)  by 


4(0  =  7T 


7  P  c‘(3)ds 


Jo  c+1 


(3.16) 


C+1  ''  4 

and  a_  €  (  —  1,0),  €  (0,1)  are  arbitrary  constants.  The  functions  ju±  are  normalized 

versions  of  Q{f^)/£  with  a  different  normalization  used  in  each  of  the  intervals  (—1,0)  and 
(0, 1).  It  is  ea.sy  to  check  that  they  are  independent  of  the  specific  choice  of  a±. 

Our  next  task  is  to  rewrite  c^  in  a  more  suggestive  form.  Given  f  =  (/i,/2,/3)  €  IR'^ 
with  fj  >  0,  the  roots  of  the  quadratic 


(/l+/2+/j)«"-(/.-/2)e-/s=0 


are  real  and  given  by 


MAf)  = 


1 


2(/i  +  /2  +  /a)  L 
1 

2(/i  +  /z  +  h) 


/l  -/2  ±  \/(/l  -  f2y  +  4/3  (/i  +  /2  +  /a) 


/i  —  /z  ±  y/ifi  +  /z  +  2/3)2  —  4/1/2 


(3.17) 


where  Ai  corresponds  to  the  minus  sign  and  A2  to  the  plus.  They  enjoy  the  properties 


-1  <  — <  A,(/)  <  0  <  AjC/)  <  <  1 . 


/i  +  /2  +  /a 


/i  +  /a  +  /a 


(3.18) 


Using  (3.17),  c®  may  be  expressed  in  the  form 


,.^J1 _ « 

e+i  i-«  < 

=  (ft  +  /I  +  «)  (« -  ^■(/'))  (« - 
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Remark.  It  is  instructive  to  compare  the  functions  with  the  eigenvalues  Xi,2{f) 

of  the  limit  fluid  equations  (1.5).  The  latter  is  a  pair  of  conservation  laws  for  /  =  (/j ,  ^2) 
of  the  type 

|a(/)  +  |b(/)  =  0.  (3.20) 

Such  a  system  is  strictly  hyperbolic  provided  that  VA(/)  is  nonsingular  and  the  roots  of 
the  characteristic  polynomial  det  (VB(f)  —  XVA(f))  =  0  are  real  and  distinct.  For  the 
case  of  (1.5),  in  the  range  of  interest  f\>0  and  /2  >  0,  the  functions  A{f)  and  B{f)  are 
smooth.  An  easy  computation  shows  that 

detV/t(/)  =  -^|=(/, +/2+v/7^)#0  (3.21) 

V/1/2 

* 

while  the  eigenvalues  are  real  and  given  by 


Note  that  Xi,2{f)  are  determined  from  a  two-vector  (/i,  /2)  while  Ai,2(/)  depend  on  a 
three- vector  (/i,  /2,  /a)  and  thus  are  not  directly  related.  However,  the  two  functions 
coincide  along  Maxwellians,  what  suggests  that  Ai,2(/^)  trace  the  eigenvalues  of  (1.5) 
along  the  fluid  dynamic  limit  £  — »  0. 

Behavior  near  the  singular  points 

Next  we  analyze  the  limiting  behavior  of  /'  in  the  neighborhood  of  the  singular  points. 
In  the  sequel  C  will  stand  for  a  generic  constant  independent  of  £. 
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Lemma  3.1.  There  are  constants  Ai±  ,  Aji  with  —1  <  Ai_  <  A14.  <  0  <  A2-  <  A2+  <  1 


depending  on  rrij  ,  Mj  but  not  £  such  that 

e  VXi_+l/ 

m3 

I «/•(«) I ,  Ti^r" 


fiitr 

TTIn 

£  ( 

e  VI-A2+/ 


for  (  e  (-l,Ai_) 
for  ^  e  (Ai  +  ,0) 
for  (  e  (0,  A2-) 
for  (  e  (A2+,1) 


As  a  consequence 


(3.23) 


~  Qif^)  0  pointwise  for  ^  £  :=  [—1,  Ai_)  U  (An-,  A2-)  U  (A2+,  1]  (3.24) 

and  uniformly  on  compact  subsets  of  £. 

Proof.  First  we  use  (2.4)  to  uniformly  estimate  the  behavior  of  c“  in  the  neighborhood  of 
the  singular  points. 

Claim:  There  exist  constants  —  1  <  Ai_  <  Ai+  <  0  <  A2-  <  A2+  <  1  such  that 


2(s  +  l) 


<  c'(.) 


c'(.)  < 


for  —  1  <  s  <  Ai_ 


for  An-  <  s  <  0 


"^3  ^  \ 


for  0  <  s  <  A2- 


c'(s)  <  - 


2(1 -s) 


for  A2+  <  s  <  1 


We  show  the  first  two  inequalities.  It  follows  from  (3.19)  and  (2.4)  that  for  s  6  (  —  1,  —h 


c'(s)  = 


/f(^)  _  /2(g)  _  /K^) 

s  + 1  3/2  1/2 

mi  M2  M3  mi 


>  ^  =•  ‘  _  A 

-  s  +  1  3/2  1/2  '  s  +  l 
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If  we  set  Ai_  =  min{  — 1/2,  —1  +  (mi/2A)},  we  see  that 


c\s)  >  -A> 


mi 


on  (-l,Ai_). 


5  +  1  2  (s  +  1) 

which  completes  the  first  inequality.  Similarly,  for  s  £  [— 2?^) 

Choose  now  A14.  =  max{— 1/2,  —mz/(2B)}  and  observe  that 

c^{s)<—  +  B<~  on  (An.,0). 
s  2s 

By  construction  the  thresholds  Ai_  ,  Ai+  have  the  property  — 1  <  Ai_  <  Ai+  <  0.  This 
proves  the  first  part  of  the  claim;  the  second  part  is  shown  in  a  similar  fashion. 

Next  we  set  a_  =  Ai_  in  (3.7)  and  use  (3.4),  (2.4)  and  the  first  inequality  in  the  claim 
to  obtain  for  — 1  <  ^  <  Ai_ 

,  QifiO) 


|^|<£exp{-irc«W..) 

c  c  ^ 


exp{-^‘ 


mi 


=  £(-i±i_) 


rn  I 

2« 


2(£  +  1) 


ds] 


Similarly,  using  (3.7)  with  o_  =  A14.  and  the  second  inequality  in  the  claim,  we  deduce 
for  Ai+  <  ^  <  0 


|«m|  <  ^  exp 


{i  f  c‘(s)ds} 


£  Ai- 
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Hence,  the  first  two  assertions  in  (3.23)  are  proved.  The  proof  of  the  other  two  is  similar 
and  is  omitted.  ■ 


§4.  Structure  of  solutions  II 

Consider  the  family  {/‘^}£>o  of  extended  solutions  to  the  problem  {Vg)  corresponding 
to  positive,  Maxwellian  data  /^.  The  functions  /'  take  the  values  f~  on  (— oo,  — 1]  and 
on  [l,oo),  they  satisfy  the  equations  (2.1  —  2.3)  for  a,  ft  6  IR.  and  enjoy  the  properties 
described  in  the  preceding  sections.  Our  objective  is  to  study  the  structure  of  the  limit 
points  of  the  family  {/^}£>o- 

Helly’s  theorem  and  the  uniform  estimates  (2.4  —  2.5)  imply  that  there  exists  a  subse¬ 
quence  (denoted  by)  {/*}  and  a  function  /  such  that  /'—»•/  pointwise  on  the  reals.  The 
limit  /  =  (/i ,  /2 1  /s )  is  a  function  of  bounded  variation,  its  components  fj  are  strictly 
positive,  and  it  satisfies  f  =  f~  on  (— oo,— 1],  /  =  /+  on  [l,oo).  Henceforth,  attention  is 
restricted  to  the  subsequence  {f^},  and  the  notation  f‘  f  will  mean  that  a  subsequence 
of  the  original  family  converges  to  the  function  /;  the  choice  of  subsequences  and  the  sense 
of  convergence  is  specified  in  context. 

The  first  lemma  establishes  that  the  sequence  {  ^  }  is  uniformly  bounded  in  .  This 

estimate  is  a  consequence  of  the  variation  bounds  amd  stronger  than  the  corresponding 
estimate  provided  by  the  balance  of  entropy  identity  (c./.  [ST]).  It  measures  hew  fast  {/'} 
relaxes  to  a  local  Maxwellian. 
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Lemma  4.1.  For  any  a,  6  €  [—1, 1],  a  <  b. 


where  the  set  S  of  exceptional  points  hcis  Lebesgue  measure  m(5)  =  0. 

Since  /  is  of  bounded  variation,  its  domain  can  be  decomposed  into  two  disjoint  subsets 
C  U  5,  such  that  on  C  all  three  components  of  /  are  continuous  while  on  S  at  least  one 
of  the  fj  is  discontinuous.  S  is  at  most  countable  and  the  right  and  left  limits  f{^+)  and 
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/(^— )  exist  at  any  ^  G  <S.  The  functions  fj  inherit  the  monotonicity  properties  of  /J;  as  a 
result  ^  =  0  is  the  only  possible  point  where  /(^— )  =  f((+)  ^  /(O  (but  this  is  excluded 
as  well  later).  So  S  includes  no  points  of  removable  discontinuity.  It  follows 

=  forfec, 

(4.2) 

M(+)  =  {SM+)M(+))''\  M(-)  =  (/,((-) foT(es, 

and  5^  C  <5  is  at  most  countable.  It  is  in  this  jjrecise  sense  that  /  is  a  local  Maxwellian.  We 
show  next  that  at  points  of  discontinuity  the  Rankine-Hugoniot  conditions  are  satisfied. 

Lemma  4.2.  For  any  point  ^  €  S 

[  /i  +  /2  +  4  \fTi~h  ]  1^^  )  +  [  “  /2  ]  =  0 

(4.3) 

-^  (  [  /l  -  /2  ]  )  +  [  /l  +  /2  ]  I 0 

Proof.  Fix  C  €  5  and  pick  any  a,  6  €  [—1,1]  —  ^  with  a  <  ^  <  b.  In  view  of  (2.1  —  2.3) 
the  sequence  {/*}  satisfies  the  identities 

-b  [/f(6)  +  mb)  +  4/1(6)]  +  a  [/f(a)  +  /|(a)  +  4/|(a)] 

+  Iftib)  -  mb)]  -  [rM  -  /!(«)]  +  f\mr)  +  mr)  +  ^mT))dr  =  O 

J  a 

-i  [f;w  -  fiw  ]  + » [  /f  («)  -  /!(<■)  ]  +  [  fU’>) + fiw  ] 

-  [/■'(») +  /l(»)]  +  /  (/.'W-ZIW)*  =  0. 

First  we  pass  to  the  limit  e  0  using  the  property  =  >7/1/2  on  C,  and  then  let  the 
points  a  — ^  (f-  and  6  ->  ^+.  In  the  combined  limit  we  obtain  the  jump  conditions  (4.3).  | 
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Since  /  may  have  discontinuities,  the  appropriate  framework  for  passing  to  the  s  0 
limit  in  (3.1  —  3.3)  is  that  of  measures.  Consider  the  functions 


Qinr)) 


dr 


(4.4) 


and  note  that  takes  constant  values  outside  [—1,1].  By  Lemma  4.1,  the  sequence 
is  uniformly  bounded  and  of  uniformly  bounded  variation  on  the  reals.  Helly’s  theorem 
implies  the  existence  of  a  subsequence,  denoted  again  by  and  a  function  of  bounded 

variation  $  such  that  $  pointwise  on  IR.  Henceforth  attention  is  restricted  to  this 

convergent  subsequence. 

Let  (p  €  C'c(IR)  be  any  continuous  function  with  compact  support  in  the  reals  and 
consider  the  Riemann-Stieltjes  integrals 


(4.5) 


Helly’s  convergence  theorem  (Natansou  [N,  VIL7])  implies  that 


< 


f 


By  the  Riesz  representation  theorem  v  may  be  viewed  eis  finite  (signed)  Borel-Stieltjes 
measures.  They  are  both  supported  on  [—1, 1],  u‘  is  generated  by  while  u  is  generated 
by  $,  the  right  continuous  version  of  $  defined  by  $(x)  =  $(x+)  (c./.  Folland  [F,  Ch.  3.5, 
Ch.  7]).  Equation  (4.6)  states  that  i/*  — *  i/  in  the  weak-*  topology  of  measures.  It  allows 
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to  pass  to  the  limit  e  ^  0  in  (3.1  —  3.3)  and  obtain 


f  h(l.V>)'di 


=  <u,lf> 


=  <u,ip> 
=  -\ 


(4.7) 


for  any  if  €  0^(11). 

The  next  lemma  characterizes  the  support  of  i/,  snppi/,  as  precisely  the  set  of  points 
where  /  is  not  a  constant  state. 


Lemma  4.3.  Let  u'  and  i/  be  as  in  (4.5).  Then 


(i)  There  are  constants  -1  <  Ai_  <  Ai^  <  0  <  X2-  <  A2+  <  1  such  that 

/_  on  (—00,  Ai_) 

f=  i  f(0)  on(Ai+,A2-) 

^  f+  on(A2+,oo) 


(4.8) 


and  supp  1/  C  [Ai _ ,  Ai+]  U  [A2- ,  A2+]. 

(ii)  (  ^  supp  u  if  and  only  if  there  exists  an  open  interval  I  containing  ^  such  that  /  is 
constant  on  I. 

(in)  S  C  supp  I/. 


Proof.  Let  Ai± ,  A2±  be  as  in  Lemma  3.1.  Integration  of  (3.3)  over  (0,^),  0  <  ^  <  A2_, 
and  use  of  (3.23)  give 


1/3(0  -  /3  (0)1 


21 


Passing  to  the  limit  £  0,  we  see  that  /j  =  /3(0)  for  €  [0,A2-).  Similar  arguments 

establish  all  other  statements  in  (4.8). 

Next  we  turn  to  the  properties  related  to  the  support  of  i/.  Recall  that  ^  suppi/  if 
there  exists  an  open  interval  /  3  ^  such  that  <  >=  0  for  any  ip  £  Cc(/).  Let  ip  be 

continuous  with  suppt^  C  (— oo,  Ai_)U  (Ai+,  A2-)  U(A2+,  oo).  Then  (4.6)  and  (3.23)  imply 

<  /'*,  p  >=  J  p  ^  —y  0  =<  p  >  as  £  —>  0. 

Hence,  suppv  C  [Ai_,  Ai+]  U  [A2-,  A2+]. 

It  suffices  to  show  (ii)  for  points  if  €  [Ai_,Ai4.]  U  [A2-,A2+].  It  follows  as  a  direct 
implication  of  the  statement: 

« 

Claim.  Let  I  be  an  open  interval  such  that  I  C  (0, 1)  or  /  C  (  —  1,0).  Then  f  is  constant 
on  /  if  and  only  if  <  u,p  >=  0  for  any  p  €  C'il)- 

One  direction  in  the  claim  is  clear  from  (4.7).  To  show  the  converse,  suppose  that  /  C  (0.1) 
(for  concreteness)  and  that  <  u^p  >=  0  for  any  p  £  Cl{I).  The  third  identity  in  (4.7) 
then  implies 

/a  t/’'  df  =  0  for  any  V’  €  Cl{I) . 

Fix  X  €  Cc{I)  with  ffX^^  =  L  Given  <f>  £  Cc{I)  define 

V’  =  J  {<f>~x{j^<i>dT))ds 

Then  0  €  Cl{I)  and  using  that  as  a  test  function  we  obtain 

^  /s  (<^  -  X  (^  <i>  dr))  =  {h  -  j ih  x)  ds)  df  =  0  for  any  €  Cc(/) . 
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It  follows  that  /s  =  Jj{f3  oil  I.  Since  is  monotone  in  (0, 1),  /s  is  constant  on 

I.  In  a  similar  fashion  it  follows  that  /i  and  /2  are  also  constant  on  I. 

This  completes  the  proof  of  part  (ii).  Part  (iii)  is  a  consequence  of  part  (ii).  g 

Along  the  convergent  sequence  f‘  f  the  functions  c'  converge  pointwise  (for  if  ^ 
±1,0)  to  a  function  c.  In  view  of  (3.19)  and  (3.8),  we  write 


1 


(/f + /I + /DU -A, (/-'))  {(-\2{n) 


1 


^  •  ~  ±  /2  ±  h)  (f  -  ^i(/))  -  ^2(f)) 

/i 


_ f2_  h 

^±1  1-f  ^ 


pointwise  on  (  —  1,0)  U  (0. 1), 


(4.9) 


and  introduce  the  function 


<7(0  = 


IL 


^€(-1,0) 
^  €(0,1) 


(4.10) 


with  a_  €  (  —  1,0),  €  (0,1)  fixed  constants. 

The  following  proposition  is  the  main  ingrecUent  of  the  analysis.  It  characterizes  points 
in  suppi^  as  being  points  of  global  maxima  (in  appropriate  intervals)  for  the  function  g. 
This  property  incorporates  admissibility  restrictions  induced  by  the  relaxation  process. 


Proposition  4.4.  Let  f  €  suppi/. 

(i)  If  f  €  suppi/  n  (0, 1)  then  g(Q  <  <jr(f)  for  any  f  €  (0, 1) 

(ii)  If  f  €  suppi/  n  (-1,0)  then  g{0  <  <7(f)  for  any  f  €  (-1,0). 


Proof.  We  show  (i).  The  proof  proceeds  in  two  steps. 
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Step  1.  rix  some  a:+  €  (0,1)  and  consider  the  functions  g‘^  defined  in  (3.8).  The 
dominated  convergence  theorem,  (4.9)  and  (2.4)  imply  that 

g‘  =  j  c‘{s)ds  g  =  I  c(s)ds  pointwise  for  G  (0, 1) . 

J  a-|.  J 

By  virtue  of  (2.4),  (3.19)  and  (3.17),  on  any  compact  [a,  6]  C  (0,1)  the  sequence  {g'^}  is 
uniformly  bounded  and  equicontinuous.  The  Ascoli-Arzela  theorem  implies  that  a  subse¬ 
quence  {<7*  }  converges  uniformly.  Since  the  whole  sequence  converges  pointwise  to  g, 


J 


c(s)  ds  uniformly  on  any  [a,  6]  C  (0, 1) 


(4.11) 


The  limit  g  is  uniformly  continuous  on  [a,  6]  C  (0,1). 


Step  2.  The  second  step  is  to  show: 

If  ^  €  (0, 1)  is  such  that  for  some  a  >  0  the  set 


.A  =  {<  e  (0, 1)  :  g{0-M  =  j  c{s)d3>a] 


(4.12) 


has  positive  Lebesgue  measme  m(A)  >  0,  then  ^  ^  supp  i/. 


To  show  (4.12),  let  ^  €  (0, 1)  be  fixed  and  let  a  >  0  and  A  be  as  above  with  Tn{A)  >  0. 
Using  (3.16)  and  (3.8),  we  write 


fl  1 


(4.13) 


Our  goal  is  to  estimate  the  integral  in  the  denominator. 
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First,  note  that  the  estimates  for  the  function  c'  in  Lemma  3.1  can  be  carried  to  the 


limit  £  -+  0  to  show:  There  are  constants  0  <  A2-  <  A2+  <  1  such  that 


m3 

2s 


<  c(s) 


c{s)  <  - 


m2 


2(l-s) 


for  0  <  s  <  A2- 
for  A2+  <  s  <  1 


Hence,  for  0  <C<A2-, 


f 

A2  — 


/  \  J  1  C 

c{s)(ls  <  In  - - >■  —00  ,  as  C  — >  0+  , 

2  '^2—' 


while,  for  A2+  <  C  <  0, 


/  c(s)  ds  <  ~  In  -  ^ - >  -00  ,  as  C  -V  1-  . 

1  ~  A2+ 

It  follows  that  points  C  near  the  endpoints  of  (0,1)  violate  the  inequality  in 
that  A  is  contained  in  some  compact  interval  [c,  d]  C  (0, 1). 

Fix  now  >  0  such  that 

-  6,  ^  implies  |  g{0  -  g{^')  |  <  ^ 

D 

Because  of  the  uniform  convergence  on  compact  subsets  (4.11),  we  can  choose 
that 

^'  €  (^  -  ^,  ^  +  ^) ,  e  <  £0  ,  C  e  imply  g‘{0  -  g‘{(')  >  ^ 
In  turn,  (4.13)  yields  for  G  (^  —  ^,  ^  +  ^) 


0  <  )  < 


< 


Ia  c+t  £ !(' 


,-(o/2r) 


(4.12),  and 


£0  >  0  such 


(4.14) 
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Let  be  a  continuous  function  with  suppv?  C  C  +  ^)-  Then  (4.5),  (4.6),  (3.14) 

and  (4.14)  together  imply 


< 


<3(/'({')) 


=  (f2(0)-f2^)  f  t^‘AO<pind^'  -  0 

JU-6,i+S) 


Thus  (  ^  suppi/  and  the  proof  of  the  statement  (4.12)  is  complete. 

It  follows  that  for  ^  €  suppi/  it  is  rn(,4)  =  0  for  any  positive  a.  Since  g  is  continuous, 
that  implies 

fl(C)  -  <7(0  =  J  ds  <  «  for  all  C  €  (0, 1) ,  a  >  0  . 

Hence,  g{Q  <  g{0  for  all  C  €  (0, 1).  | 

The  maximization  properties  of  g  at  points  in  the  support  of  u  provide  information  on 
the  structure  of  the  limit  function  /.  In  particular,  a  weak  form  of  the  Lax  shock  conditions 
is  satisfied  at  any  point  of  jump  discontinuity. 

Corollary  4.5.  Let  0  €  [—1, 1],  ^  <  0- 

(a)  If  n  supp  u  then 

<  =  Ai(/(0)  for^<0, 

(4.15) 

e  =  A2(/(0)  fore>0. 

(b)  If  ^  €  <S  then  /  satisfies  at  ^  the  jump  conditions  (4.3)  and  the  inequalities 

A,(/(^+))  <  ^  <  A.(/(^-))  fore<0, 

(4.16) 

A2(/(^+))  <  i  <  Hm-))  forC  >0. 
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(c)  If  €  suppi/n(-l,0)  then  Ai(/(^+))  =  Xiif {('-))  =  and  for  6  €  (6^'): 


e  =  \,{m) 

if  0eC 

(4.17) 

Ai(/(0+))  =  ^  =  AU/(M) 

if  6  €S 

(4.18) 

If  ('  e  suppr/  n  (0,1)  then  A2(/(^+))  =  \2ifiC-))  =  i'  and  for  d  €  (tC)- 


^  =  Wm  a  dec 


(4.19) 


A2(/(^+))  =  0  =  \2ifi0-))  if  ees 


(4.20) 


Proof.  The  proof  is  presented  for  ^  €  (0, 1);  the  case  of  (  —  1,0)  is  .similar. 

« 

Consider  the  function  g  in  (4.10)  restricted  to  the  interval  (0,1).  Note  that  g  is  contin¬ 
uous  and  g{Q  —00  as  (  — »  O-f  or  ^  >  1—.  Since  c  is  of  bounded  variation 

lim  — £^11=  lim  ^  ■  f  c{s)  ds  —  c{(±) .  (4.21) 

C-«±  C  -  ^  C  -  ^  A 

Thus  ^  exists  and  is  continuous  at  points  of  C,  while  only  the  right  and  left  derivatives 
exist  at  points  of  S.  Let  ^  €  supp  u  fl  (0, 1).  By  Proposition  4.4  the  function  g  achieves  its 
global  maximum  on  (0, 1)  at  ^ 

giO  <  g{0  for  any  C  €  (0, 1) . 

First  we  show  (a)  and  (b):  At  ^ 

c(e+)<o,  c(c-)>o. 
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from  where,  using  (4.9),  (3.18)  and  noting  that  ^  >  0,  we  conclude 

e  -  Hm+))  >  0 ,  e  -  A2(/(e-))  <  0 .  (4.22) 

If  ^  6  C  then  (4.22)  leads  to  (4.15).  If  ^  6  <S  then  the  jump  conditions  (4.3)  are  satisfied 
at  by  Lemma  4.2,  and  (4.22)  gives  (4.16). 

Next  we  show  (c).  The  points  6  supp  i/f!  (0, 1),  ^  are  both  global  maxima  for 
g.  Therefore  g(^)  =  g{(')-  We  claim: 

9{^)  =  <7(0  for  any  9  6  (OC')  •  (4.23) 

If  (4.23)  is  violated  at  some  point,  there  exist  a,  b  with  ^  <  a  <  6  <  0  such  that 

gia)  =  g(b)  =  g{^) ,  <7(0  <  <7(0  for  a  <  0  <  6. 

At  the  points  a,  b  we  have 

t 

A2(/(a+))  <a<  A2(/(a-)) 

A2(/(6+))  <  b  <  A2(/(6-)) 

On  the  other  hand,  it  follows  from  Lemma  4.3  and  Proposition  4.4  that  /  is  constant  in 

9 

the  interval  (a,  6)  and  thus  A2(/(a+))  =  A2(/(6— )).  The  inequedities  then  imply  b  <  a, 
which  contradicts  a  <  b-,  hence  (4.23)  follows. 


28 


This  establishes  that  g  stays  constant  at  its  maximum  value  on  [^,  Therefore  c(^+)  = 
c(^'— )  =  0  and  c{d—)  =  c(0+)  =  0  for  ^  ^  In  turn,  these  yield 

0  =  A2(/(^))  for  0€C 

X2{f{e+))  =  d  =  X2ifid-))  for  0es 

^  =  X2ifU+)) 

e'  =  A2(/(e'-)) 

and  the  proof  is  complete,  g 

Behavior  of  solutions  in  the  wave  fans 

It  follows  that  the  region  where  /  is  noncoiistant  consists  of  two  disjoint  closed  intervals; 
fxi  =  C  (—1,0)  associated  with  the  first  characteristic  speed  Ai(/)  of  (1.5),  and 

=  [02,  h]  C  (0, 1)  associated  with  the  second  characteristic  speed  A2(/).  Each  of  Ix^  or 
1x2  could  be  empty  or  consist  of  just  a  single  point.  The  function  /  takes  constant  values 
on  the  complement  of  Ixi  U  /a^,  while  it  has  the  behavior  indicated  in  Corollary  4.5  at 
points  of  Ja,  or  Ixj- 

Our  next  objective  is  to  obtain  a  fuller  description  of  the  behavior  of  /  on  the  two 
wave  fans.  The  following  proposition  does  not  use  the  specific  form  of  (2.7),  but  only  its 
structural  properties  and  relations  (4.17  —  4.20).  It  is  thus  convenient  to  work  with  the 
abstract  form  of  (2.7) 


29 


and  the  associated  weah  form 


-bAifib))  +  aAifia))  +  Bifib))  -  Bifia))  + 


A(/(s))  ds  =  0 


a,beTR.  (4.24) 


Consider  the  strictly  hyperbolic  system  (3.20)  with  A(f),  B{f)  at  least  twice  continu¬ 
ously  differentiable  and  VA(f)  nonsingular.  Let  A,(/)  be  its  eigenvalues,  li{f)  the  left 
eigenvectors  and  ri(/)  the  right  eigenvectors.  They  are  connected  through  the  relations 

V5(/)r.(/)  =  A,(/)V^(/)r.(/) 


(4.25) 


/.(/)  •  VB(/)  =  A.(/)  hif)  •  VAif) 


and  satisfy  the  normalization  conditions 


li(f]-VA(f)r,{f)=S,,.  (4.26) 

Proposition  4.6  Suppose  that  =  [a/fc,6i]  is  a  full  interval,  Uk  <  bk- 

(i)  For  each  ^  €  [ajt,6jt)  such  that  VAfc(/(^-|-))  •  rk{f{(+))  7^  0  it  is 

fc-o”i>o  h  ~  "  VAfc(/(e+))  •  rfc(/(e+)) 

(ii)  For  each  ^  €  (ajt,^*]  such  that  VAjt(/(^-))  •  rjfc(/(^-))  /  0  it  is 

» J”i<0 1 
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Proof.  We  show  (i).  Fix  ^  €  [a*:,  6fc)  and  let  h>  0  such  that  ^  +  /i  €  /a*-  The  weak  form 
(4.24)  taken  between  the  points  and  ^  +  h—  gives 


[  -e  VA(/(a))  +  VB(/($+))  ]  (/({  +  h-)  -  /({+))  (4.29) 

=(«+'■)  [A(m + h-)  -  A{m+)  -  vAm+))  (/(? + »>-)  -  /(?+))  ] 

-  [B(/(«  +  h-)  -  B(/(?+)  -  VB{f{i+))  (/(e  +  /.-)  -  /({+))  ] 

r^+h 

The  increment  (/(^  +  h—)  —  /(if+))  is  expanded  in  the  basis  of  right  eigenvectors 


^{h)  :=  /(^  +  h-)  -  /(^+)  =  5;]u;.(/0n(/(^+)  (4.30) 

X 

Note  that  for  a  function  of  bounded  variation  a;(/i)  ->  0  as  h  — >  O-h,  and  that  by  (4.26) 


>^i{h)  =  um+)-'^A{m+))u{h). 


(4.31) 


Taking  the  inner  product  of  (4.29)  with  li{f{^+)  and  using  (4.25),  (4.31)  and  the  Taylor 
expansion,  we  obtain 


1  +  A, (/(«+))  I  u,i(/.)  =  0(Mfc)l")  +  0(  f  M^)]  ds  )+0{k  |u.{MI)  .  (4.32) 

On  account  of  Corollary  4.5  and  the  strict  hyperbolicity  property  of  (1.5),  the  coefficient 
(— ^  +  Ai(/(^+))  is  nonzero  for  i  /  /;  but  vanishes  for  i  =  k. 

Next,  using  (4.17  —  4.20)  and  the  Taylor  expansion  of  A;t)  we  see  that 


Afc(/(e  +  h-))  -  \km+)}  =  h  (4.33) 

=  vA,(/(e+))  •  im  +  ft-)  -  /(«+))  +  o(MA)p) . 
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If  we  set  gk  =  VA(t(/(^4-))  •  rk{f(^+)),  Qk  ^  0  by  hypothesis,  and  use  (4.33),  (4.30)  and 
relations  (4.32)  for  i  fc  we  arrive  at  the  estimate 


gk^k{h)  -  fe  =  0(5]k(MI)  +  0(KA)p) 


=  O(MA)P)  +  0(  /  ds)+Oih  |a;(A)|)  ,  (4.34) 


Adding  (4.32)  for  i  ^  k  with  (4.34)  gives 


^{h)  ■  =  \gkOJk{h)  -  /il  +  ^  |wi(/i)| 


=  o(  (Mfc)i  +  h)  |^(/oi )  +  0(  r  Ms)l  ds ) 

=  0(  (|a,(A)|  +  hMk]  )+0(j  9(^)  ds )  +  0{h^) 


(4.35) 


Since  ui{h)  — >  0  as  /i  — >  0+,  we  can  choose  6  sufficiently  small  so  that  for  0  <  h  <  6 


^(h)<Ch^+C  /  ^is)ds 
Jo 


The  integral  inequality,  in  turn,  yields 


0<ip{h)<C'h^  {ov0<h<6 


and  thus 

u}i{h)  n  r  ■  /  1  T  ^ 

hm  — ; —  =  0  for  i  ^  k ,  hm  — ; —  =  —  . 

fc— 0+  h  h— 0+  h  gk 

This  shows  (4.27).  Part  (ii)  is  proved  similarly.  | 

Proposition  4.6  implies  that  /  has  right  and  left  derivatives  at  any  point  ^  which  is  not 
an  accumulation  point  of  S.  If  such  a  point  ^  belongs  to  C  then  /  is  Lipshitz  there,  and 
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if,  in  addition,  it  is  an  interior  point  of  /a*  then  /  is  differentiable  there.  It  also  completes 
the  picture  regarding  the  structure  of  the  wave  fans.  There  are  the  following  cases: 

(i)  If  Ja*  consists  of  a  single  point  then  the  solution  is  a  shock  wave  satisfying  the  weak 
form  of  the  Lax  shock  conditions  (4.16). 

(ii)  If  is  a  full  interval  of  points  in  C  the  solution  is  a  k-rarefaction  wave  (provided  that 

•  T’fc  0  on  /aj,  which  is  anyway  necessary  for  rarefactions). 

(iii)  In  general  /a*  consists  of  an  alternating  sequence  of  shock  waves  and  k-rarefaction 
waves  such  that  each  shock  adjacent  to  a  rarefaction  from  one  side  is  a  contact  discontinuity 
on  that  side. 

The  emerging  picture  is  that  typical  for  small  data  of  strictly  hyperbolic  but  not  genuinely 
nonlinear  systems.  It  will  be  further  simplified  in  Section  6,  by  using  the  geometry  of  the 
shock  curves  for  (1.5). 

§5.  Self-similar  limits  and  shock  profiles 

Our  next  task  is  to  discuss  the  relation  between  self-similair  fluid  dynamic  limits  and 
shock  profiles  for  the  Broadwell  model.  To  this  end  fix  ^  a  point  of  discontinuity  for  /, 
and  note  that  /(^— )  ^  /(C+)  are  both  M2ixwellian  states  and  satisfy  (4.3). 

Consider  a  sequence  of  points  {.^e}  with  the  property  —>■  to  be  specified  later. 
Define  functions  Vj(Q  =  fji^e  +  eC)  the  new  variable  — oo  <  C  <  oo.  This  accounts 
for  a  shift  of  the  original  solution  and  the  introduction  of  the  stretched  variable  C  The 
uniform  estimates  (2.4  —  2.5)  imply 
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(5.1) 


0<m,  </;(^,+eC)<  Mj, 

TV^f‘{^.+e-)  =  TV^f]{-)<  K,. 

Using  Helly’s  theorem  and  a  diagonal  argument  we  establish  the  existence  of  a  subsequence 
and  a  function  v  =  (vi,V2,V3)  such  that 

fjUs  +  ^C)  ^j(C)  pointwise  for  —  oo  <  C  <  oo  .  (5.2) 


The  sequence  defined  in  (4.4)  is  of  uniformly  bounded  variation.  Upon  restricting  to 
further  subsequences,  if  necessary,  we  also  have 


$ 


«  QU'M) 


dr  —*  $(0  pointwise  for  —  oo  <  ^  <  oo  ,  (5.3) 


-  $ 


0  . 


The  following  lemma  shows  that  the  phase  shift  {^4}  can  be  arranged  so  that  v  is  a 
shock  profile  for  the  Broadwell  system  connecting  the  states  /(^— )  with  /(^+). 

Proposition  5.1  There  exists  a  choice  of  the  sequence  {^4}  such  that  v(C),  defined  by 
(5.2),  is  continuously  differentiable  and  satisfies  the  differential  equations 

dui  dvi 


dv2  dv2  . 

1  .  . 


(5.4) 


and  the  boundary  conditions 


u(-oo)  = /(^-),  v{+oo)  =  f{^+)  . 


(5.5) 
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Proof.  The  proof  hcis  two  steps.  The  first  step,  to  show  (5.4),  is  independent  of  the 
particular  choice  of  the  sequence 

We  evaluate  (2.1)  between  the  points  -h  £(  and  ^  and  obtaui 


-(G  +  --  C)  /i'(6  +  £  C)  +  ^  /i(^)  +  /f  (6  +  e  0  -  AW 

r^«+^C  2 


(5.6) 


+  /  Air)dT=-  Q{nr))dT 

Je  ^  Je 


Observe  that,  along  a  subsequence  and  ais  £  0, 


/  A{T)dT  /  Air)  dr 
Je  Je 


while 


1  y?<+eC 

-  /  Q(r(r))dr 
£  Je 


1  1  r*'^*'* 

=  -/  Q{nT))dT  +  -  Qinr))dT 

£  Je  £  J?. 

=  ($^(^,)  -  $'(0))  +  Q{r{^.+es))ds 

Jo 

^  ($°  -  $(0))  +  [  Q{v{s))  ds . 

Jo 


Passing  to  the  limit  in  (5.6),  we  deduce  that  for  0,  C  €  IR 


-«»i  (C)  +  f .  (C)  -  Q(v{s))ds  =  -9/, (9)  +  /,  (9)  +  /  /,  (r)  dr  +  (*”  -  4(9)) .  (5.7) 


In  a  similar  fashion  from  (2.2  —  2.3)  we  obtain  the  relations 


-iMO  -  t-2(C)  -  Q(v(s))ds  =  -9/2(9)  -  /2(9)  +  f  MT)dT  +  ($“  -  *(9)) ,  (5.8) 
-(vM)+lJ^^Q{v{s))ds  =  -«M0)  +  I  Mr) dr- (5.9) 
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Since  the  above  equations  hold  for  any  C  G  IR-i  both  the  left  and  the  right  hand  sides  are 
constant.  That  is,  there  are  six  equations  hidden  in  (5.7  —  5.9);  three  for  /,  and  three  for 
V  describing  the  internal  structure  of  the  shock  profile.  The  functions  Vj  are  continuously 
differentiable  and  satisfy  the  system  of  ordinary  differential  equations  (5.4). 

We  proceed  to  show  (5.5).  In  preparation,  note  that  (5.7  —  5.9)  may  be  combined  to 
obtain 

-((v,{0  +  2v,{0)  +  v^i0  =  A 

(  V2{a  +  2  V,i0 )  -  V2(0  =  B  (5.10) 

-^^(C)  =  -^(^3^(C)-t^i(C)r^2(C)), 

% 

where  the  constants  ^1,  B  are  calculated  from  the  right  hand  sides  of  (5.7  — 5.9),  as  0  » if±, 
and  are  given  by  the  formulas 

^  =  -«(/,«+) +2 /,«+)) +/,«+)  =  -e  (/■«-)+ 2 /,«-))  +  /,(«-) 

B  =  -?  (/2K+)  +  2 /,«+))  -  /,«+)  =  (/2(e-)  +  2/3(?-))  -  /,(?-) 

Also  recall  from  (4.2) 

/3''(«+)  -/l«+)/2«+)  =  0,  h\(-)  -  Mi-)  hd-)  =  0. 

The  equilibria  (1x1,027^3)  of  (5.10)  are  determined  by  solving  the  algebraic  system 

— (uj  +  2u3)  +  ui  =  a 

— ^  ( XX2  +  2  U3  )  —  U2  — (5T1) 
U3^  —  Ui  U2  =  0 
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Substituting  the  first  two  expressions  in  the  third,  we  see  that  U3  is  determined  from  the 
roots  of  the  quadratic 

(l-^")tX3"+(^  +  2^U3)(5  +  2^t/3) 

and  then  ui  ,  U2  are  obtained  by  solving  the  linear  equations.  Therefore  (5.11)  has  at  most 
two  real  solutions,  and  as  a  result  the  equilibria  of  (5.10)  are  precisely  /(^— )  and  /(^+). 

We  turn  now  to  (5.5).  To  fix  ideas,  suppose  that  ^  E  (0, 1)  and  that  fj  is  monotone 
increasing  (for  some  j)  in  (0, 1)  but  may  be  increasing  or  decreasing  in  (  —  1.0).  Let  d  >  0 
be  fixed.  Given  any  bounded  Interval  J  3  0  and  provided  that  £  is  chosen  sufficiently 
small,  the  function  f  C)  is  monotone  increasing  on  J  and 

/;(^-^)</;(6+eC)</;(^  +  ^),  for  cej. 

Passing  to  the  limit  ^  >  0,  we  deduce  that  Vj  is  monotone  increasing  on  any  J,  and  thus 
also  on  (— cxDjOo),  and  that 

<  Vj(-oo)  <  vj(0)  <  vj(+oo)  <  fj(^+) .  (5.12) 

Both  the  inequalities  and  the  monotonicity  properties  are  reversed  when  we  start  with  a 
component  fJ  that  is  monotone  decreasing  in  (0, 1).  In  either  case  the  limits  v(— 00)  and 
t;(+oo)  exist  and  are  finite. 

For  ^  €  5  it  is  /3(^-)  ^  /3(^+);  otherwise  /(^-)  =  /(^+)  by  the  analysis  of  (5.11). 
Suppose  for  concreteness  that  <  /3(^+)  and  fix  a  state  1/3(0)  such  that 

hU-)<V3{0)<hU  +  ). 
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Since  /I  is  monotone  near  ^  and  /I  /a  pointwise,  for  each  sufficiently  small  e  we  can 
choose  such  that  flUe)  =  ^3(0).  The  sequence  {i^^}  has  the  property  ifi  —*■  With 
this  choice  of  define  the  cissociated  v  by  (5.2).  Then  v  satisfies  (5.10)  and  the  limits 


u(— oo)  and  i;(+oo)  exist  and  satisfy  (5.12).  Since  f{^—)  and  f{^+)  are  the  only  equilibria 
of  (5.10)  and  v(0),  by  virtue  of  selection,  is  not  an  equilibrium  state,  we  conclude  that 

/(^-)  ^  u(-oo)  ^  i;(+oo)  =  /(^+) 


and  the  proof  is  complete.  | 

§6.  The  limit  fluid  equations 

« 

In  this  section  we  use  the  special  properties  of  the  limit  fluid  equations  (1.5)  to  complete 
the  description  of  tlie  wave  fans.  The  properties  used  are  the  genuine  nonlinearity  of  the 
characteristic  fields  of  (1.5)  and  the  geometry  of  the  shock  curves.  Although  they  are  both 
discussed  in  Caflisch  [C],  we  give  an  independent  presentation  for  completeness  and  to 
account  for  extended  differences  in  notation. 

The  eigenvalues  Ai(/),  \2if)  satisfy  the  characteristic  polynomial 

ifi  +  f2  +  V /i/z)  —  ifi  —  f2)  —  y/ fih  =  0  (6.1) 


and  are  given  by  (3.22)  or  (3.17)  for  /j  =  y/JTh-  The  corresponding  right  eigenvectors  are 
ezisily  computed 
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Ai  + 1 
Ai-1 
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i  =  l,2. 


We  differentiate  (6.1)  in  order  to  express  in  terms  of  A^,  and  use  the  resulting 

expressions  to  compute 

VA  .  r  -  ~  {  (/l  +  f2+  4\//i/2)  Ai  -  (/i  -  /2)  \ 

2\//i/2  \  (/l  +  /2  +  \/fif2)  2\i  —  {fi  —  f2)) 

The  denominator  in  (6.2)  is  positive  for  the  positive  eigenvalue  A2  and  negative  for  A}.  On 
the  other  hand,  (3.17)  yields 


(/l  -  /2)  -  1/1  -  /2|  ^  /l  -  /2 
2  (/i  +  f2  +  vZ/i/a )  /i  +  /2  +  4\//]/2 
(/l  -  /2)  +  1/1  -  /2I  ^  /l  -  /2 
2  (/i  +  /2  +  s/fTH)  f\  +  f2  +  ^y/Tih 


It  follows  from  (6.2)  and  (3.18)  that  VA,-  •  7-^  >  0  for  i  =  1,  2. 

Since  (1.5)  is  genuinely  nonlinear,  the  pos.sibility  of  contact  discontinuities  is  ruled  out 
at  least  for  weak  shocks.  For  the  case  of  strong  shocks  we  need  to  study  the  geometry  of 
the  shock  curves. 

The  shock  curves  are  defined  by  solving  the  Rankine-Hugoniot  conditions 


—s  (  ui  +  2  U3  )  +  Ui  =  —s  (  +  2  U3  )  +  Vj 


-S  (  U2  +  2  ll3  )  -  U2  =  -5  (  U2  +  2  U3  )  -  772 


(6.3) 


U3^  -  Ui  U2  =  U3^  -  V2 


=  0 


for  the  states  u  =  (ui,U2,U3)  and  v  =  {vi,V2,V3)  and  the  shock  speed  s.  To  this  end  fix 
one  state,  say  v,  with  V3  =  y/viV2  and  consider  the  increments  [u^]  =  Ui  —  Vi.  Then 

2s  (U3]  =  (1  -5)[ui]  =  -(1  +  s)[u2]  (6.4) 

[W-l]  ([us]  +  2  773)  =  [ui]  [772]  +  V2  [u,]  -f  771  [7x2]  (6.5) 
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Substituting  (6.4)  into  (6.5)  yields 


(1  +  3s^)[u3]  -  2  [(ui  +U2  +U3)-S^  -  (ui  -  U2)5  -  U3] 


=  0 


If  [U3]  =  0  there  is  no  shock.  Thus,  using  (6.1),  we  obtziin  a  representation  of  the  shock 
curve  parametrized  by  the  speed  s  in  the  form  of 


2(ui  +  U2  +  yJV\V2) 
1  +  3s2 


(3-\,{v))  (S-\2(V))  ■ 


(6.6) 


together  with  (6.4).  As  a  consequence,  if  the  shock  is  a  contact  discontinuity  (at  either 
side)  then  the  strength  of  the  shock  [u]  =  0. 

This  observation  excludes  the  possibility  of  contact  discontinuities  and  simplifies  the 
structure  of  the  solutions  for  (1.5)  considerably.  We  conclude  then  that  each  wave  fan  is 
either  a  single  rarefaction  or  a  single  shock  which  satisfies  the  Lax  shock  conditions. 
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